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ABSTRACT 

We report the identification of 61.45 (711.2 p-Hz) oscillations, with amplitudes 

of 62.6-pmag, in KIC 4768731 (HD 225914) using Kepler photometry. This relatively 
bright (V=9.17) chemically peculiar star with spectral type A5 Vp SrCr(Eu) has pre¬ 
viously been found to exhibit rotational modulation with a period of 5.21 d. Fourier 
analysis reveals a simple dipole pulsator with an amplitude that has remained stable 
over a 4-yr time span, but with a frequency that is variable. Analysis of high-resolution 
spectra yields stellar parameters of = 8100 ± 200 K, logg = 4.0 ± 0.2, [Fe/H] = 
-1-0.31 ±0.24 and usint = 14.8 ± 1.6 kms“^. Line profile variations caused by rotation 
are also evident. Lines of Sr, Cr, Eu, Mg and Si are strongest when the star is brightest, 
while Y and Ba vary in anti-phase with the other elements. The abundances of rare 
earth elements are only modestly enhanced compared to other roAp stars of similar 
Teff and logg. Radial velocities in the literature suggest a significant change over the 
past 30 yr, but the radial velocities presented here show no significant change over a 
period of 4yr. 

Key words: asteroseismology - stars: abundances - stars: chemically peculiar - 
stars: individual (KIC 4768731) - stars: magnetic field - stars: oscillations - techniques: 
photometric 


1 INTRODUCTION 


The chemically peculiar Ap stars are a spectroscopic 
subclass of A-type stars which exhibit strong enhance- 
men ts of one or m ore elements, notably Cr, Eu, Si and 
Sr (iMorganI 1 19331). Strong m agnetic fields in Ap stars 
were found by Babcoclj lll947t) and their presence appears 
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ubiquitous with field stren gths reaching several kilogauss 
dDonati fc Landstreed 120091 ) . As a group the Ap stars ro¬ 
tate more slowly than n ormal A-type stars due to magnetic 
breaking (IStepieii|l200Cll ) and exhibit photometric and spec¬ 
tral variability due to inhomogeneous surface distributions 
of element al abundances concentrat ed into spots by the mag¬ 
netic field llWolflllig's^ : ISmithlll996l ). 

The rapidly oscillat ing Ap (roA p) stars are a subset of 
the magnetic Ap stars (lKurtdll982l ). These exhibit short- 
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Table 1. Catalogue of roAp stars. The first two columns are the HD number and other name if available. The effective temperature 
(K) is obtained from the Balmer line profiles when available. The luminosity, T/Lq is obtained from Tgff and the surface gravity when 
available or from the radius. The p rojected rotational velocity, usini (kms“^) is given. The averaged magnetic field, B (kG), from 
iBychkov, Bychkova & Madeil ll2009l) are shown; values in italics are other measures of B. The rotation period, Prot is determined from 
the light curve. A maximum of four pulsation frequencies is given. 
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timescale variations with periods between 5 and 25 minutes 
and amplitndes up to 0.01 mag. The roAp stars are rela¬ 
tively rare with 61 known to d ate (Table fill, compared t o 
around 2000 known Ap stars iRenson fc ManfroidI l2009ll . 
Most have been fonnd using ground-based photometry of 
already known Ap stars. Some, which do not exhibit pho¬ 
tometric variations detectable from the gronnd, have been 


fonnd by spect roscopic studies looking for radial velocit y 
variations (e.g. iKochukhov et aDl2002l : lElkin et all 1200,^ 1. 
Eleven of the 61 known roAp stars have now been fonnd 
using the S uper WASP archive with E-band amplit udes 
>0.5mmag (|Holdsworth et ^l2014al : iHoldsworthllioi^ l. 

The launch of the Kepler spacecraft enabled the search 
for photometric oscillations with amplitudes well below 
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those detectable from the ground. Seven stars previously 
classified as Ap stars were observed by Kepler to search 
for rapid oscillat ions, but only one w as found to exhibit 
roAp pulsations (iBalona et al. l2011bl i. Two further previ¬ 
ously unknown Ap stars were identifi ed from their pulsa¬ 
tions and subsequ ent spectral analyses llBalona et aLllioilal : 
iKurtz et al.ll201ll '). In the sp ectroscopic survey of 11 7 A and 
F stars observed by Kepler, iNiemczura et ah (l2015l i identi¬ 
fied one new Ap star (KIC 476873 1), and re-classified two of 
the stars by iBalona et al.l (l2011bl) as not Ap (KIC 8750029 
and KIC9147002). All the Kepler-discovered roAp stars 
have oscillation amplitudes below O.lmmag, while one ini¬ 
tially discovered by SuperWASP located in the Kepler fiel d 
has an amplitude of over 1 mmag (iHoldsworth et ani2014bl) . 
The star discussed here, KIC 4768731, is brighter and has a 
lower pulsation frequency than the other known roAp stars 
in the Kepler field. 


2 KIC 4768731 

Prior to the Kepler mission, KIC 4768731 (HD 225914, 
BD -1-39 3919) w as a rather anonymous V = 9.17 sta r 
dHog et al.ll2000l ) with a spectral type A7 (ICanno ^[liii). 
The star is part of the double system WDS 19484-1-3952 with 
a 12th magnitude star (KIC 4768748) at 12" dMason et al.l 
[20 o 3). There has been very little change in separation or 
position angle over a period of 74 yr, suggesting that the 
two stars might be a com mon proper motion pair. Recent 
prop e r motion catalogues (iRoeser, Demleitner fc Schilbachl 
120101 : IZacharias et al.[ 12013 ) are inconclusive in this re¬ 
gard, with both stars having proper motions that agree to 
within the error bars. Un fortunately, the Hipparcos parallax 
dPerrvman fc ESj^ll997h of tt = 8.68 ± 7.80 mas is too un¬ 
certain to provide any meaningful distance constraint. The 
distance of ~350 pc estimated from the stellar parameters 
(see Sect. [Tjl would imply a physical separation of at least 
4000 au. However, the Kepler Input Catalog estimates Tefs 
= 4220 K and logg = 1.55 for the companion star, suggest¬ 
ing that this it a background K-type giant star. In Sect. [6] 
spectroscopic observations confirm this. 


3 SPECTRAL CLASSIFICATION 

KIC 4768731 was included in the sample of several hun¬ 
dred A- and F-type stars to be observed at high-resolution 
as part of the Kepler Asteroseismic Science Consortium 
(KASC ) spectroscopic follow- up of stars with Kepler light 
curves dNiemczura et aHl2015h . A spectrum of this star was 
taken on 2011 July 13 using the fibre-fed High Efficiency 
and Resolution Me rcator Echelle Spectrograph (HERMES; 
iRaskin et al.ll20lll) . Fig.[T] presents the spectrum at classi¬ 
fication resolution. This spectrum is a smoothed version of 
the high-resolution spectrum downgraded to a resolution of 
1.8A to match that of the spectral standard^. 

The H 7 line and the strength of the metal lines in gen¬ 
eral match that of the A5V standard, but for an A5V star the 


^ Spectra of the MK standards were 

obtained at R.O. Gray’s website: 

http: //stellar .phys. appstate. edu/Staiidards/stdl_8 .html 



Wavelength (A) 

Figure 1. Spectrum of KIC 4768731 compared to the A5 V and 
A7V MK standards. Note the prominent Sr II lines at 4077A and 
4216A and the abnormal Call K line. The spectra have been offset 
for clarity. 


Ca II K line is shallow and abnormally broad. This abnormal 
shape is due to abundance stratification within the star’s 
atmosphere (see Sect. 17311 . The Cai 4226A line, however, 
does not look abnormal. Sril lines are particularly strong 
(see the 4077 and 4216A lines). Cr lines are also strong: the 
Cm 4111A line that is blended in the redward core-wing 
boundary of HJ is pronounced, and this is confirmed with 
the Cm 3866A line. The Cm 4172A line is mixed in with 
the Fe-Ti 4172-9A blend, but the 4172A component is no¬ 
ticeably stronger. 

Other features to look for in the spectrum of Ap stars 
include strong Eu and Si lines (for an atlas of Ap star spec¬ 
tra, see lCrav fc Corballvll200'9l ). The Eull 4205A line is en- 
hanced, but the Eull 4130A line is blended with the Sill 
doublet. That doublet of Sill (4128 and 413lA) is strong 
and, unlike in the comparison spectra, is more of a blended 
feature than a separated doublet. This is not a product 
of rotation, as KIC 4768731 is clearly not a rapid rotator, 
but could reflect a strong Eu line. The Sin 3856A line ap¬ 
pears strong, but other Sill lines (4002, 4028 and 4076A) 
are not pronounced. These findings lead to the classification 
A5Vp SrCr(Eu). 


4 KEPLER OBSERVATIONS 

IPebosscher et aP ll201ll) classified KIC 4768731 as exhibit¬ 
ing rotational modulation with a frequency of 0.191206 d“^ 
(5.23 d) and an amplitude of 5.407 mmag. The AAVSO Vari¬ 
able Star Index fVSX: IWatson|[2006h noted this star as ACV: 
(a suspected CVn variable) with a period of 5.21 d and 
magnitude range of 10 mmag based on Kepler light curves. 
The VSX entry by R. Jansen was dated 2011 May based on 
a phase plot of Kepler data from JD 2454964 to 2454998, 
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Figure 2. Kepler pixel map for the SC data overlaid on the 
2MASS J-band image. The dark grey area contains the pixels 
used in the standard PDCSAP extraction and in light grey are all 
the pixels downloaded from the spacecraft. The areas highlighted 
in dashed and dotted lines are those used here to extract light 
curves of KIC 4768731 and KIC 4768748, respectively. 


noting also that the rise time from minimum to maximum 
light is 50 per cent of the variable’s period. 

One month of Short Cadence (SC) data was acquired by 
KASC for this star in Q2.1 (2009 June 20 - 2009 July 20), 
but with a 2-d gap due to a safe mode event. The star was, 
nevertheless, observed nearly continuously in Long Cadence 
(LC) mode throughout the Kepler mission from QO to Q17 
(2009 May 02 - 2013 May 2013). A periodogram of the SC 
data shows the previously reported rotational modulation, 
plus a weak higher frequency signal at 61.45 d“^. The vari¬ 
ability of KIC 4768731 obtained from the Kepler light curves 
is discussed in detail in Sect. [5] 

The presence of the nearby companion star 
(KIC 4768748) could potentially affect the Kepler light 
curve of KIC 4768731. Fig. [2] shows the SC pixel map 
overlaid with the 2MASS J-band image. The standard 
PDCSAP light curve extraction does not include the compan- 
ion star, but it is ver y close. Therefore, using pyke tools 
dStill fc Barclay! I 2 OI 2 I '). we extracted two light curves; one 
containing just the pixels centred on KIC 4768731 and the 
other centred on the companion star KIC 4768748. The light 
curve of KIC 4768748 contains neither the rotational mod¬ 
ulation nor the high-frequency signal, while the light curve 
of KIC 4768731 show the previously reported rotational 
modulation plus the weak higher frequency signal. 
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Figure 3. Top panel: An amplitude spectrum for the Q1 to Q17 
Kepler long cadence data out to nearly the Nyquist frequency 
where only the low-frequency rotational variations can be seen, as 
a consequence of the scale. Bottom panel: An amplitude spectrum 
after a high-pass filter has removed the low-frequency variation 
showing the Nyquist alias of the 61.45-d~^ pulsation frequency 
triplet. Note the three orders of magnitude change in the ordinate 
scale from mmag in the top panel to pmag in the bottom panel. 


5 FREQUENCY ANALYSIS OF THE KEPLER 

DATA 

KIC 4768731 was observed by Kepler in all quarters QO-17 
in LC, and during one month, Q2.1, in SC. We use these data 
to study the pulsation frequencies, amplitudes and phases 
for the full 4-yr data set. There is no Nyquist limitation with 
Kepler LC data, although there is a reduction in amplitude 
as a consequence of the pulsat ion periods being shorter than 
the i ntegration times. See iMurohv. Shibahashi fc Kurt j 
(l2013tl for an explanation. Fig.O shows an amplitude spec¬ 
trum of the QO-17 data where only the low-frequency peaks 
generated by the rotational variation can be seen in the top 
panel. We ran a high-pass filter that removed the rotational 
variations from the light curves, and the lower panel of Fig.[3] 
shows the amplitude spectrum of those filtered data. There 
is a clear frequency triplet centred at 12.51 d“^. A closer look 
at these peaks shows them to be Nyquist aliases of higher 
frequency peaks near 61.5 d“^, as is evident from the SC 
data. 

The rotation frequency found by fitting the highest peak 
in Fig. [ 3 ] to the LC data by nonlinear least-squares gives 
j/rot = 0.1919622 ± 0.0000003 d■^ or P^ot = 5.209358 ± 
0.000009 d. 

Fig.|4]shows an amplitude spectrum of the Kepler Q2.1 
SC data in the top left panel, where a clear peak is seen at 
61.45 d~^. After pre-whitening by the highest peak, the ro¬ 
tational sidelobes at — r'rot and i^iJ-r'rot become evident on 
the lower left panel. This frequency triplet is equally spaced 
and is evidence of an oblique dipole pulsation mode. This 
shows that the 12.51 d“^ peak seen in Fig.[3]is a Nyquist alias 
of this peak, since 61.45 - z/lc = 12.51, where z/lc = 48.94 
is the Kepler LC sampling frequency. The top right panel 
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Figure 4. Top left panel: An amplitude spectrum for the Q2.1 SC showing the pulsation frequency near 61.45 d“^. Bottom left panel: 
The SC data after pre-whitening by where the rotational sidelobes are evident. Top right panel: The same frequency range for the 
QO—17 LC data where the oblique dipole pulsation triplet can be seen. Bottom right panel: The LC data after pre-whitening by ui, where 
the rotational sidelobes are evident. The residual amplitude around is the result of frequency or phase variability of the pulsation 
frequency over the 4-yr time span of the data set, as is seen in Fig. [5] 


of Fig. 2] therefore shows the same section of the ampli¬ 
tude spectrum for the entire QO-17 Kepler LC data set, 
where the frequency resolution is the highest available. Of 
course, the amplitude is reduced because the integration 
times of the LC data are longer than the pulsation period. 
From iMur iy lioTi) we calculate the reduction factor to 
be AjAo = sm{n/n)/{n/n) — 5.5, where n is the number of 
data points per oscillation cycle; this explains the lower am¬ 
plitude in the LC data when compared to the SC data. The 
bottom right panel then shows the amplitude spectrum of 
the LC data after prewhitening by z^i, where the rotational 
sidelobes can be seen, along with residual amplitude around 
u\. This residual amplitude is a consequence of frequency 
variability over the 4-yr data set, as we show below. 

This triplet seen in Fig. U for both the SC and LC data 
is eqnally split by exactly the rotation freqnency to better 
than Icr. This is the signature of an oblique dipole pulsation 
mode. The three frequencies have been fitted to the data by 
a combination of linear and nonlinear least-squares. Table[2] 
shows the results where the frequency splitting has been 
set to be exactly the rotation frequency. A time zero point 
was selected to set the phases of the sidelobes equal, and 
we found that the phase of the central frequency is also 
equal, showing that the triplet is the result of pure amplitude 
modulation. This is typical of pure oblique dipole pulsation. 
We then find that the time of pulsation maximum, when all 
three frequencies of the triplet have the same phase, to = 
2455697.97056, coincides with maximum brightness of the 
rotational light modulation. That is shown by the phase of 
the rotational frequency, near —7t, with that time zero point. 
Because we have fitted a cosine function, and the data are 
in magnitudes, the phase at maximum rotational brightness 
is 7t radians. This can also be seen in the light curve, which 
we do not show here. 

Maximum pulsation amplitude at rotational light mini¬ 


mum is expected for dipole oblique pulsation where the spots 
that produce the rotational light variations are aligned with 
the magnetic and pulsation poles. For KIC 4768731 the re¬ 
lationship differs by n radians from this simplest case, sug¬ 
gesting that the spots, magnetic field pole and pulsation 
pole lie in the same plane, but are not completely aligned. 
This is consistent wit h the improved oblique pulsator model 
dBigot fc Kurt j[201ll i. 

We can constrain t he pulsation geometry within the 
oblique pulsator model dKurt j1 19821 : iBigot fc Kurt j I 2 OIII I 
with the amplitudes of the components of the frequency 
triplet: 

tan i tan [3 = (A+i -|- A-\)/Aq, 

where i is the rotational inclination, p is the inclination 
of the magnetic pole to the rotation axis, A+i, A_i and 
Ao are the amplitudes of the rotational sidelobes and cen¬ 
tral frequency. From this equation we obtain tan i tan P = 
0.226 ± 0.031 and tan i tan/I = 0.233 ± 0.016 from LC and 
SC data, respectively. Using the stellar radius and usini 
obtained from spectral analysis (Sect. [T]) and the rotation 
period, we obtain values of i = 401 ^q ° and P = 151,0 °- 
The angle between the magnetic and rotation axes is there¬ 
fore relatively small. In addition, i + P < 90° and only one 
pole is seen, hence the simple, almost sinusoidal rotational 
light variations. 

To study the residual amplitude around seen for 
the LC data, we fitted vi to 20-d sections of the data by 
least-squares and examined the time series of amplitude and 
phase. Fig.[3shows that the amplitude is stable over the 4-yr 
time span, and the phase is variable. Since phase and fre¬ 
quency variability are indistinguishable, we interpret this as 
frequency variability. This is the source of the residual am¬ 
plitude in the amplitude spectrum that is unresolved from 
!zi as seen in the bottom right panel of Fig.|4l Such frequency 
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Table 2. A least-squares fit of the rotation frequency and the pulsational dipole triplet for KIC 4768731 to the full 4-yr Kepler LC 
data set. The zero point of the time scale is BJD2455697.97056, which coincides with pulsation amplitude maximum and rotational light 
maximum. This is typical for an oblique pulsator. 


labels 

frequency 

d-i 

amplitude 

|i.mag 

phase 

radians 

^'rot 

rotation 

0.1919622 ±0.0000003 5368.3 ± 3.8 - 

-3.0211 ± 0.0007 


pulsation - long cadence (LC) data 


^rot 

FI - 1 - Frot 

61.255837 ±0.000088 
61.447800 ±0.000008 
61.639762 ±0.000054 

0.8 ±0.2 

9.3 ±0.2 

1.3 ±0.2 

1.112 ±0.222 
1.117 ±0.020 
1.112 ±0.147 


pulsation - short cadence (SC) data 


1^1 l^TOt 

Fl + Frot 

61.255837 ±0.000088 
61.447800 ±0.000008 
61.639762 ±0.000054 

6.5 ±0.7 

62.2 ±0.7 

8.0 ±0.7 

2.842 ±0.115 

3.057 ±0.012 
2.842 ± 0.092 



0 300 600 900 1200 

BJD-2454963 (d) 

Figure 5. Top panel: pulsation amplitude of ui = 61.4478 d~^ as 
a function of time for 20-d segments of the QO-17 KIC 4768731 
data. The amplitude is stable, but may have a secular decrease 
over the 4 yr. Bottom panel: pulsation phase for the same seg¬ 
ments. There is clear variability in the pulsation phase. There is 
an apparent oscillation with a period of about half the Kepler or¬ 
bital period, i.e. 186 d, which may be instrumental in origin. The 
larger, longer-term variations are most probably astrophysical in 
origin. 

va riability is known in many other roAp stars, as is discussed 
bv iHoldsworth et al.l ll2014bfl for the roAp star KIC 7582608. 
The cause of such frequency variability is not yet known. 


6 SPECTROSCOPIC OBSERVATIONS 

To examine the abundances in detail and to investigate line 
profile variations as the star rotates, we obtained 33 high- 
resolution spectra. 

Two high-resolution {R=85 000) spectra covering the 
wavelength range 3780-9000A were obtained using the fibre- 
fed High Efficiency and Resolution Mercator Echelle Spec¬ 


trograph (HERMES; I Raskin et al.l [20llh mounted on the 
1.2-m Mercator Telescope at the Roque de los Muchachos 
Observatory, La Palma. The data were reduced using the 
HERMES DR^ pipeline software. 

A further two high-resolution spectra with a resolu¬ 
tion R=67 000 covering the wavelength region 3660-7360A 
were obtained with the 2.5-m Nordic Optical Telescope 
(NOT) using the Fibre-fed Echelle Spectrograph (FIES; 
iTelting et al.[[20l3l . These spectra were extracted with the 
bespoke data reduction package, FIESTOOI0. 

A series of 26 high-resolution spectra were obtained in 
2013 August using the coude-echelle spectrograph attached 
to the 2.0-m telescope at the Thiiringer Landessternwarte 
(TLS) Tautenburg, Germany. The spectrograph has a re¬ 
solving power of i? = 63 000 and exposures cover the wave¬ 
length range from 4720-7360A. The spectra were reduced 
using standard ESO MIDAS packages. All spectra were cor¬ 
rected in wavelength for individual instrumental shifts by 
using a large number of telluric O 2 lines. 

Table [ 3 ] gives the dates of the observations, along with 
the rotational phase {4>iot) at the time of the midpoint of 
the exposures and measured radial velocity (see Sect. ED)- 
The rotational phase is obtained from the fractional part of 
(RJH —2455697.97056) x 0.1919622, where phase zero is the 
time of maximum light. 

A TLS spectrum of the companion star, KIC 4768748, 
obtained on 2014 June 6, revealed that the star is indeed 
a late-K giant. In addition, the radial velocity of this star 
is —2.9 ± 0.1 kms“^, compared to —12.2 ± 0.2 kms“^ for 
KIC 4768731 (See Sect. 17.211 . Thus, we conclude that the 
two stars are not physically associated, KIC 4768748 is a 
background giant. 


7 SPECTROSCOPIC ANALYSIS 

A detailed spectroscopi c analysis of the 2011 HE RMES spec¬ 
trum was presented in iNiemczura et all (|2015h who identi- 


^ http://www.mercator.iac.es/instruments/hermes/hermesdrs.php 
® http://www.not.lac.es/instruments/fies/fiestool/ 
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Table 3. Log of spectral observations of KIC 4768731. Date and BJD are given for the mid-point of the 
exposure. 


Instrument 

Date 


Exp Time (s) 

S/N 

BJD (TDB) 

0rot 

RV (kms-i) 

HERMES 

2011 Jul 13 22:53:56 

1900 

87 

2455756.45765 

0.227 

-12.0 

EIES 

2013 Aug 04 

00:43:06 

2000 

102 

2456508.53370 

0.597 

-12.3 

EIES 

2013 Aug 06 

00:19:54 

1748 

112 

2456510.51759 

0.978 

-11.8 

HERMES 

2013 Aug 10 

03:26:10 

1800 

80 

2456514.64693 

0.771 

-12.2 

TLS 

2013 Aug 14 

22:59:41 

1800 

80 

2456519.46185 

0.695 

-12.3 

TLS 

2013 Aug 14 

23:30:42 

1800 

73 

2456519.48339 

0.699 

-12.3 

TLS 

2013 Aug 15 

00:01:43 

1800 

70 

2456519.50493 

0.703 

-12.4 

TLS 

2013 Aug 15 

00:32:44 

1800 

63 

2456519.52646 

0.708 

-12.4 

TLS 

2013 Aug 15 

01:41:49 

1800 

58 

2456519.57444 

0.717 

-12.3 

TLS 

2013 Aug 15 

02:12:50 

1800 

49 

2456519.59598 

0.721 

-12.3 

TLS 

2013 Aug 15 

02:43:51 

1800 

46 

2456519.61752 

0.725 

-11.7 

TLS 

2013 Aug 15 

21:54:18 

2400 

108 

2456520.41643 

0.879 

-12.1 

TLS 

2013 Aug 16 

01:58:32 

2400 

85 

2456520.58604 

0.911 

-12.7 

TLS 

2013 Aug 16 

21:44:06 

2400 

95 

2456521.40934 

0.069 

-12.4 

TLS 

2013 Aug 17 

01:59:27 

2400 

71 

2456521.58667 

0.103 

-12.4 

TLS 

2013 Aug 17 

22:59:50 

2400 

89 

2456522.46192 

0.271 

-12.3 

TLS 

2013 Aug 21 

19:54:11 

2400 

76 

2456526.33295 

0.014 

-11.9 

TLS 

2013 Aug 22 

20:10:07 

2400 

72 

2456527.34400 

0.208 

-11.9 

TLS 

2013 Aug 23 

19:51:48 

2400 

76 

2456528.33127 

0.398 

-12.2 

TLS 

2013 Aug 24 

00:50:49 

2417 

70 

2456528.53891 

0.438 

-12.2 

TLS 

2013 Aug 26 

20:58:21 

2400 

67 

2456531.37743 

0.983 

-12.4 

TLS 

2013 Aug 26 

22:25:29 

2400 

78 

2456531.43794 

0.994 

-12.4 

TLS 

2013 Aug 26 

23:06:30 

2400 

73 

2456531.46642 

1.000 

-12.3 

TLS 

2013 Aug 27 

01:42:35 

2400 

60 

2456531.57481 

0.021 

-12.3 

TLS 

2013 Aug 27 

20:42:41 

2400 

45 

2456532.36653 

0.172 

-12.1 

TLS 

2013 Aug 27 

21:24:35 

2400 

98 

2456532.39563 

0.178 

-12.2 

TLS 

2013 Aug 27 

22:07:30 

2400 

95 

2456532.42543 

0.184 

-12.3 

TLS 

2013 Aug 27 

22:48:31 

2400 

84 

2456532.45392 

0.189 

-12.3 

TLS 

2013 Aug 27 

23:29:32 

2400 

90 

2456532.48240 

0.195 

-12.4 

TLS 

2013 Aug 28 

00:10:33 

2400 

76 

2456532.51088 

0.200 

-12.3 

TLS 

2014 Jun 17 

23:10:15 

2400 

86 

2456826.46545 

0.628 

-12.4 

TLS 

2015 May 10 22:55:57 

2400 

84 

2457153.45487 

0.398 

-12.2 

TLS 

2015 May 11 01:07:39 

2400 

75 

2457153.54733 

0.416 

-12.1 


Notes: The rotational phase (</>rot) is given by the fractional part of {BJD — 2455697.97056) X 0.1919622, 
where phase zero is the time of maxim um light. The signal-to-noise (S/N) ratio obtained using the 
DER_SNR algorithm llStoehr et alJl^OSli . 


Table 4. Basic stellar parameters of KIC 4768731. 


Parameter 

Value 

Units 

Feff 

8100 ± 200 

K 

logs 

4.0 ± 0.2 


[Fe/H] 

+0.31 ± 0.24 


^'^icro 

0.5 ± 0.3 

kms“^ 

i/sini 

14.8 ± 1.6 

kms“^ 

M 

2.11 ± 0.27 

M© 

R 

2.39 ± 0.68 

R© 

logL 

1.34 ± 0.25 

L© 

Mv 

1.43 ± 0.67 

mag. 

Sp. Type 

A5Vp SrCr(Eu) 



Note: Mass (M) and rad i us (R ) are estimated using the 
Torres. Andersen fc Gimenej ll201Clli calibration. 


7.1 Magnetic field 


The ratio of the strengths of the Fell 6147.7A and 6149.2A 
lines can be used to estimate the mean magnetic field 
modulus (H) (iMathvs fc Lanj|l992l) . From the individual 
spectra we obtain a mean value for the relative intensih- 
cation of the Fell 6147.7A line with respect to the Fell 
6149. 2A line, AW\/W\, of 0.08 ± 0.03, where AW\ = 
Wa(6147.7) — VFa (6149.2). This yields a mean magnetic field 
modulus, (H ), of 2.7 ± 0.8 kG usin g the empirical relation¬ 
ship given in lMathvs fc Lan3 (Il992l ). The error bar includes 
both the scatter in the observed values of AW\/W\ and the 
uncertainty in the calibration for {H). The individual values 
for AW\/W\ show no significant variation with rotational 
phase. 


For the magnetic modulus obtained above the split¬ 
ting of th e 6149.2A Zeeman doublet components would be 
-^O.IS A l|Mathvs et ahlflOO^ . However, there is no sign of 
fied KIC 4768731 as an Ap CrSrEu star. Table [4] summarizes any splitting in our spectra due to the relatively large v sin i 

the basic stellar parameters obtained for this star. of this star. 


© 0000 RAS, MNRAS 000, 000-000 






















Barry Smalley et al 


- 11.4 

- 11.6 
- 11.8 
f" -12 

I -12-2 

^ - 12.4 

- 12.6 
- 12.8 
-13 


Figure 6. Radial velocity measurements show no trend with ro¬ 
tational phase. The filled circles indicate the values obtained in 
2011, 2014 and 2015 and show that they are not significantly 
different from those from 2013 (open circles). 

7.2 Radial velocity variations? 

A radial velocity (RV) for each spectrum was determined 
by cross-correlation with a synthetic spectrum covering the 
wavelength range 5000-5800A using the parameters and 
abundances obtained from the detailed spectroscopic anal¬ 
ysis. The barycentric radial velocities are given in Table [3] 
and the formal uncertainties are 0.2 kms“^. The average 
radial velocity from the spectra obtained in 2013 is — 12.2 
± 0.2 kms“^. However, iFehrenbach fc Burnagd lll990l) re¬ 
ported a value of — 29±1.2 kms“^ based on four observa¬ 
tions obtained on 1984 June 27, 1986 August 2, 30 and 31. 
They used the Fehrenbach Objective Prisms on the Schmidt 
telescope at the Observatorie de Haute-Provence and the ex¬ 
ternal error is given as 3.5 kms~^. The difference between 
the two epochs of ~17 kms~^ suggests that the star could 
be a spectroscopic binary. The HERMES spectrum taken in 
July 2011 has a radial velocity that is consistent with those 
taken during August 2013. In June 2014 and May 2015 we 
obtained further TLS spectra and their RV values are con¬ 
sistent with the previous years (Fig. [6|. Hence, we find no 
significant change in RV over a period of 4yr. 

7.3 Spectral line profile variations 

The individual spectra were examined for evidence of spec¬ 
tral line profile variations with rotational phase. The Ha and 
H/1 Balmer lines and the NaD lines do not show any varia¬ 
tion with phase. Similarly, the HERMES and FIES spectra 
of the Ca H and K lines do not show any variation with 
rotational phase, but have the unusual line profiles associ¬ 
ated with abundance stratificat ion in the atmosphere dBabell 
1 19941 : iRvabchikova et "ni 120021 ) . The sharp cores of Ca H 
and K lines can be recreated using a step-shaped stratified 
abundance profile with [Ca/H] =4-1.0 dex for layers deeper 
than log(T5ooo) = —1.0 and [Ca/H] = —2 dex for higher lay¬ 
ers (Fig. (7)). In contrast, a global enhancement of [Ca/H] = 
4-0.5 dex is unable to reproduce the observed line profile. 

The HERMES and FIES spectra of the Srii 4077.7A 
line shows considerable variation with rotational phase 
(Fig. [8]). The line is strongest at c/rot = 0.978 and weak¬ 
est at 0rot = 0.597. This is consistent with the weakest line 



Wavelength (A) 


Figure 7. The Ca H and K lines in the FIES spectrum (grey 
line) taken on 2013 August 6 shows the distinctive pointed core 
due to element stratification. The solid line shows a synthetic line 
profile created using a step-shaped stratified abundance profile, 
while the red dashed-line shows a normal non-stratified synthesis. 



Wavelength (A) 

Figure 8. Variations in the profile shape of the Sr II 4077.7A line 
with rotational modulation phase. The solid black line was taken 
nt (prot = 0.978, the solid red line at (/>rot = 0.597. The dashed 
green and dotted blue lines were taken at (/> = 0.771 and tpiot = 
0.227, respectively. 


strength occurring at or around />rot = 0.5. The other two 
spectra are very similar, but their phases with respect to 
/>rot = 0.5 are —0.271 and 4-0.273, respectively. This sug¬ 
gests that there is a region of enhanced Sr crossing the stellar 
meridian at the time of maximum light {(prot = 0.0). 

Fig. [9] shows the variation of the Y ll 4883.7A line pro¬ 
file with rotational phase. There is a clear absorption fea¬ 
ture moving redward with increasing phase, first appearing 
on the blue edge of the line profile at around i/rot =0.1 and 
disappearing around (prot = 0.9. The feature crosses the cen¬ 
tre of the line around (prot =0.5, coinciding with the time 
of minimum light in the rotation cycle. Similar behaviour 
is seen in the Yll 4900.1 and Ban 6141.7A lines. This is in 
anti-phase to the variation of the Sril line. 

Fig. [To] shows the variation of the Eu ll 6437.OA line 
profile with rotational phase. The line is quite weak around 
(prot — 0.4. By (prot = 0.7 an absorption feature can be seen 
in the blue wing of the profile. As phase increases the feature 
moves redward across the line profile, reaching a maximum 


2011 


2015 


2014 


0 0.2 


0.4 0.6 

Phase (() ,) 


0.8 
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Table 5. Elemental abundances obtained for KIC 4768731 from TLS spectra. Abundances are given in the form 
logA(El) = log(A^Ei/A^H) + 12 obtained from spectra taken on six different nights and rotational phases (0rot)- fi 
is the maximum number of lines for an element used in the abundance averages, which can be slightly fewer at different 
rotat ional phases. L i ne-to- line standard deviation error bars are given when n > 2. In column 9 the Solar abundances 
from lAsplund et al.l ll2009l) are given in for reference. 


El 

n 

Aug 14 

(0rot = 

-15 
: 0.71) 

Aug 15 

{<pTOt — 

-16 
: 0.89) 

Aug 16 

("^rot = 

-17 
: 0.09) 

Aug 23 

(0rot = 

-24 
: 0.42) 

Aug 26 

(<^rot = 

-27 
: 1.00) 

Aug 27-28 
=0.19) 

Solar 

C 

19 

8.45 

± 

0.33 

8.26 

± 

0.27 

8.47 

± 

0.36 

8.51 

± 

0.28 

8.43 

± 

0.34 

8.46 

± 

0.37 

8.43 

N 

1 

8.16 



8.38 



8.06 



8.20 



7.94 



8.15 



7.83 

O 

5 

9.37 

± 

0.22 

9.10 

± 

0.51 

9.14 

± 

0.41 

8.62 

± 

0.55 

9.16 

± 

0.42 

9.40 

± 

0.21 

8.69 

Na 

3 

7.13 

± 

0.06 

7.08 



7.21 

± 

0.12 

6.81 

± 

0.31 

7.22 

± 

0.12 

6.77 



6.24 

Mg 

8 

7.83 

± 

0.18 

7.87 

± 

0.21 

7.82 

± 

0.20 

7.73 

± 

0.17 

7.86 

± 

0.20 

7.84 

± 

0.16 

7.60 

A1 

2 

7.47 



7.43 



7.45 



7.04 



7.34 



7.15 



6.45 

Si 

34 

7.53 

± 

0.27 

7.66 

± 

0.24 

7.75 

± 

0.21 

7.45 

± 

0.26 

7.74 

± 

0.23 

7.56 

± 

0.22 

7.51 

S 

9 

7.34 

± 

0.43 

7.56 

± 

0.49 

7.57 

± 

0.25 

7.02 

± 

0.26 

7.27 

± 

0.33 

7.22 

± 

0.34 

7.12 

Ca 

18 

6.81 

± 

0.21 

6.77 

± 

0.20 

6.88 

± 

0.21 

6.80 

± 

0.18 

6.81 

± 

0.20 

6.83 

± 

0.22 

6.34 

Sc 

7 

2.45 

± 

0.10 

2.47 

± 

0.35 

2.46 

± 

0.33 

2.66 

± 

0.45 

2.30 

± 

0.35 

2.53 

± 

0.45 

3.15 

Ti 

42 

5.08 

± 

0.29 

5.00 

± 

0.28 

5.13 

± 

0.25 

5.12 

± 

0.24 

5.17 

± 

0.24 

5.19 

± 

0.27 

4.95 

V 

8 

5.10 

± 

0.30 

5.17 

± 

0.30 

5.19 

± 

0.28 

5.14 

± 

0.30 

5.22 

± 

0.28 

5.21 

± 

0.24 

3.93 

Cr 

111 

7.36 

± 

0.18 

7.50 

± 

0.18 

7.51 

± 

0.17 

7.21 

± 

0.16 

7.55 

± 

0.17 

7.38 

± 

0.19 

5.64 

Mn 

17 

5.77 

± 

0.35 

5.82 

± 

0.36 

5.72 

± 

0.32 

5.79 

± 

0.28 

5.84 

± 

0.34 

5.69 

± 

0.30 

5.43 

Fe 

142 

7.79 

± 

0.18 

7.83 

± 

0.18 

7.83 

± 

0.17 

7.80 

± 

0.19 

7.81 

± 

0.18 

7.83 

± 

0.18 

7.50 

Co 

6 

5.78 

± 

0.25 

6.11 

± 

0.37 

6.14 

± 

0.33 

5.84 

± 

0.24 

6.28 

± 

0.26 

5.74 

± 

0.18 

4.99 

Ni 

31 

6.30 

± 

0.34 

6.37 

± 

0.36 

6.33 

± 

0.37 

6.18 

± 

0.28 

6.24 

± 

0.35 

6.39 

± 

0.31 

6.22 

Cu 

1 

3.09 



3.05 



3.12 



3.34 



3.63 



3.08 



4.19 

Zn 

1 

4.73 



4.70 



4.53 



4.68 



4.69 



4.72 



4.56 

Y 

9 

3.05 

± 

0.24 

2.74 

± 

0.39 

2.90 

± 

0.25 

3.12 

± 

0.21 

2.82 

± 

0.20 

2.79 

± 

0.35 

2.21 

Zr 

3 

3.60 



4.67 

± 

0.77 

4.36 

± 

0.69 

3.49 



3.48 



4.27 

± 

0.58 

2.58 

Ba 

2 

2.50 



1.81 



1.87 



2.52 



1.80 



2.05 



2.18 

La 

2 

1.60 



1.34 



1.91 



0.77 



1.23 



1.75 



1.10 

Nd 

6 

2.34 

± 

0.15 

2.24 

± 

0.02 

2.31 

± 

0.14 

2.19 

± 

0.24 

2.40 

± 

0.02 

2.28 

± 

0.06 

1.42 

Eu 

3 

2.65 

± 

0.31 

2.50 

± 

0.11 

2.33 

± 

0.36 

1.93 

± 

0.50 

2.47 

± 

0.34 

2.20 

± 

0.26 

0.52 



Wavelength (A) 

Figure 9. Changes in line profile strength and shape with rota¬ 
tional phase for the Y II line at 4883.7A. The rectified spectra have 
been offset so that their continuum levels correspond to their rota¬ 
tional phases. The nearby Cm line at 4884. fiA shows little shape 
variation with phase. The vertical black and grey lines indicate 
the line centres and the ±15 kms“^ usinf widths, respectively. 


around 0rot = 0.0 when it is close to the line centre. The 
feature then moves to the red wing and weakens, disappear¬ 
ing around 0rot = 0.2. Similar behaviour is seen in other Eu 
lines. The variations in the Eu lines are in anti-phase with 
those for Y and Ba lines, but in phase with the Sr line. 

To further investigate line strength variations with ro- 



Wavelength (A) 

Figure 10. Changes in line profile strength and shape with ro¬ 
tational phase for the Full line at 6437. fiA. The rectified spec¬ 
tra have been offset so that their continuum levels correspond to 
their rotational phases. The nearby Cal line at 6439. lA shows lit¬ 
tle shape variation with phase. The vertical black and grey lines 
indicate the line centres and the ±15 kms“^ usini widths, re¬ 
spectively. 

tational phase, we selected six nights in 2013 August with 
multiple TLS spectra. For each night the spectra were co¬ 
added to produce a single nightly spectrum with signal- 
to-noise (S/N) ratios in the range 100 - 150. For each of 
the six spectra abundances were obtained using the sam e 
spectral synthesis methodology as iNiemczura et al.l ll2015ll . 
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Table 6. Rare earth elements in KIC 4768731. Abundances are 
given in the form logA(El) = log(AfEi/fVH) + 12. A colon (:) 
indicates a tentative detection, where other lines of the element 
have upper limits below the given value, n is the number of lines 
used to obtain the abundance averages, except in the cases of 
upper limits or tentative detections where it is the number of lines 
search ed. In column 4 the Solar abundances from lAsplund et alj 
1 I 2 OO 9 I ') are given in for reference. 


El 

n 

logA(El) 

Solar 

LS/ 

2 

1.40 

1.10 

Ce 

6 

<2.5 

1.58 

Pr 

3 

3.2: 

0.72 

Nd 

6 

2.25 ± 0.12 

1.42 

Sm 

2 

2.2: 

0.96 

Eu 

3 

2.45 ± 0.11 

0.52 

Gd 

4 

3.0: 

1.07 

Tb 

3 

<2.5 

0.30 

Dy 

4 

<3.3 

1.10 

Er 

2 

3.5: 

0.92 

Tm 

2 

<2.5 

0.10 

Yb 

6 

<3.5 

0.84 

Lu 

5 

<2.0 

0.10 


The results of the spectral fitting are given in Table O The 
overall abundance patt ern is consistent with that found by 
iNiemczura et all (l2015l l and with the spectral classification 
presented in Sect. [3] Cr is considerably enhanced, Eu is mod¬ 
estly enhanced, while Ca and Si are relatively normal. V and 
Co are also considerably enhanced and there is general mod¬ 
est -1-0.3 dex) overabundance of most other transition- 
group elements, with the exception of Sc, which is depleted 
by over 0.5 dex. Several elements. Mg, Si, Cr, Co and Eu, 
have abundance variations in phase with rotation - they are 
strongest when the star is the brightest. In addition, the 
same behaviour was noted for the Sr line above. The Y and 
Ba lines, on the other hand, vary in anti-phase with rotation 
- they are strongest when the star is the faintest. Other ele¬ 
ments, including Ca, Ti, V, Mn and Fe, show little variation 
with phase or the results are inconclusive. 

To search for the presence of other rare earth elements, 
the individual TLS spectra from 2013 were co-added to pro¬ 
duce a single spectrum with a S/N of 245. This spectrum 
was visually searched for lines due to various rare earth el¬ 
ements. With the exception of the already detected La, Nd 
and Eu, only tentative detections or upper limits were ob¬ 
tained (see Table O. 

The relatively large, and variable, scatter in the stan¬ 
dard deviations of some of the abundances obtained is prob¬ 
ably the result of the surface inhomogeneities and vertical 
stratification present in the atmosphere of this Ap star. A 
detailed spectral analysis taking into account surface inho¬ 
mogeneities and abundance stratification will be presented 
in Niemczura, Shulyak et al. (in preparation). 


8 NON-ADIABATIC OSCILLATION 
MODELLING 

The mechanism responsible for exciting the oscillations ob¬ 
served in roAp stars is still not fully understood. The opacity 
mechanism acting on the hydrogen ionization region leads 


to the excitation of high radial-order acoustic pulsations in 
models of cool Ap stars with fully radiative envelopes, where 
convection is assumed to be supp ressed by the strong mag¬ 
netic field llBalmforth et al]l200in . This mechanism has been 
shown to provide a promising explanation for the observed 
pulsations in roAp stars with p ulsation frequencies below 
the acoustic cutoff llCunhall2002l L However, in a fraction of 
the known roAp pulsators, the observed frequencie s are too 
high to be driven in this way (ICunha et al.l |2()L^ and an 
alternative mechanism needs to be con sidered, such as the 
effect of turbulent pressure suggested bv ICunha et al.l (l2013fl 
and recently fo und to be a likely e xplanation for pulsations 
in an Am star (lAntoci et al.ll^Ol^ L 

The oscillation frequency observed in KIC 4768731 is 
below the acoustic cutoff frequency derived from models cov¬ 
ering the region of the HR diagram where the star is located. 
To check if the opacity mechanism can drive the observed 
pulsation, and further constrain the global parameters of the 
star, we carried out a linear, non-adiabatic stability analy¬ 
sis of a grid of models covering that region of the HR di¬ 
agram. Evolutionary tracks were produc ed with the MESA 
stellar evolution code JPaxton et al.llioi^ l for stellar masses 
in the range 1.6-2.4 Mq, in intervals of 0.05 Mq, and with 
initial mass fraction of hydrogen and helium oi X = 0.70 
and Y = 0.28, respectively. Mixing beyond the convective 
core during the main-sequence has been cons idered, taking 
the co nvective overshoot parameter defined in IPaxton et al.l 
(l2013ll (their equation (9)) to be /ov = 0.016. Only mod¬ 
els with Teff in the range 7850-8350 K were considered for 
the stability analysis (vertical dashed lines in Fig. 1111) . The 
mass, luminosity, and effective temperature extracted from 
the MESA models were used to generate the equilibrium mod¬ 
els necessary for the non-adiabatic computations. 

The non-adiab atic analysis followed closely that of 
ICunha et al.l ll2013ll . except that we have considered only 
their ‘polar’ models (models with fully radiative envelopes) 
as these are the models in which roAp-type pulsations are 
found to be driven by the opacity mechanism. Following 
that work, for each set of mass, luminosity, and effective 
temperature, we have considered four different case stud¬ 
ies. In the first (or standard ) case the equilibrium model 
has a surface helium abundance of Yaur = 0.01 and an 
atmosphere that extends to a minimum optical depth of 
Vmin = 3.5 X 10“®. Moreover, in the pulsation analysis a 
fully reflective boundary condition is applied at the surface 
for this case study. The other three case studies are ob¬ 
tained by swapping these properties, one at the time to: 
Yaur = 0.1, Tmin = 3.5 X 10“’’, and a transmissive boundary 
condition. Together, these four case studies cover the main 
uncertainties in the modelling. If modes with frequencies 
similar to that observed in KIC 4768731 are found unstable 
in at least one of the case studies, then the corresponding 
model parameter set is identified as unstable at the observed 
frequency. 

Models found to be unstable at the observed frequency 
are marked by filled-circles in Fig. 1111 Inspection of this 
figure shows that there are two independent regions within 
the part of the HR diagram explored in which pulsations 
with the observed frequency are predicted to be excited. 
In the lower luminosity region, the modes excited at the 
observed frequency are of low radial-order, characteristic of 
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Figure 11. HR diagram showing the results of the non-adiabatic 
oscillation modelling. The cross indicates the location and Icr un¬ 
certainty of KIC 4768731. The solid lines are the evolutionary 
tracks (starting from the ZAMS) covering masses from 1.6 Mq 
( bottom) to 2.4 (top) in steps of 0.05 Mq. The vertical dashed- 
lines indicate the range in temperatures searched. Models found 
to be unstable at the observed frequency are marked by filled- 
circles. In the lower luminosity region, the modes excited are of 
low radial-order, while in the higher luminosity region they are of 
high radial-order. 


^-Scuti stars, while in the higher luminosity region they are 
of high radial-order, characteristic of roAp stars. 

The non-adiabatic analysis performed here does not in¬ 
clude the direct effect of the magnetic field on pulsations, but 
only its indirect effect through the suppression of convection 
in the stellar envelope. However, it is known from studies 
addressing the direct effect of the magnetic field that the 
coupling of acoustic and magnetic waves in the outer layers 
of roAp stars leads to wave energy losses t hrough the dissi¬ 
patio n of acoustic waves in the atmosphere dSousa &: Cunha 
20081 ) and mag netic waves in the interior (ICunha fc GoughI 
20001) . In fact, ISaid (l2005l) has shown that energy losses 
through the latter process are particularly significant for 
low radial-order modes. These modes are then stabilised, 
explaining their absence in roAp stars. If this is the case 
for KIC 4768731, then the lower luminosity unstable region 
shown in Fig. [TT] is likely to be spurious, resulting from 
the non-inclusion of the direct effect of the magnetic field 
on pulsations in the non-adiabatic computations. Our non- 
adiabatic calculations thus point to KIC 4768731 being a 
relatively evolved and luminous star. 


9 DISCUSSION AND CONCLUSIONS 

Given the identification of KIC 4768731 as an Ap star, 
the pulsation is almost certainly of the roAp type, and 
with a frequency of 61.45 d“^ (corresponding to a period 
of 23.43 min) this star only just misses ou t on being the 
longe st period roAp star known (HD 177765; lAlentiev et al.l 
|2012|. 23.56 min). Other low-frequency roAp stars, discov¬ 
ered from their low-amplitude radial velocity va riations, in- 
clude HD 116114 (21 min) and 0 CrB (16.2 min) llElkin et al.l 
l2005bl : iKurtz. Elkin fc Mathvdl2007^ . All four of these stars 
occupy a similar location in the HR diagram, around TbH 


Figure 12. Abundance patterns for KIC 4768731 (filled circles), 
HD 177765 (open circles), HD 116114 (open squares) and /3 CrB 
(crosses). Upper limits are indicated by downward pointing ar- 
rows. Abundances a re given relative to the Solar values from 
lAsplund et al.l ll2009h . 


= 8000 K and logg = 4.0, with KIC 4768731 having the 
shortest rotation period. 

Spectral lines f rom rare e arth elements are usually 
strong in Ap stars llWolfll [l983h . However, we find a dis¬ 
tinct lack of significant overabundances in rare earth ele¬ 
ments compared to those found in oth er roAp stars of sim i- 
lar effective temperatur es: HD 177765 jAlent iev et al.ll2012l ). 
HD 116114 and /3 CrB dR.vabchikova et al.ll2004l ) fFig. 1121) . 


A similar conclusion was found for another , cooler, Kepler 
roAp star, KIC 10195926 dElkin et al.ll20l3) . 

The Mg, Si, Cr, Co, Sr, and Eu lines in KIC 4768631 
show line strength variations in phase with rotation - they 
are strongest when the star is the brightest. The Y and Ba 
lines, on the other hand, vary in anti-phase with rotation 
- they are strongest when the star is the faintest. Other 
elements, including Ca and Fe, show little variation with ro¬ 
tational phase or the results are inconclusive. The pulsation 
analysis has shown that the pulsation axis is in the same 
plane as the abundance spots that lead to rotational light 
maximum, but not completely aligned. The pulsations are 
coincident with light maximum, hence with Sr ll maximum, 
and out of phase with Y ll maximum. Where the pulsation 
axis, magnetic poles and spots are located is problematic, 
but they appear to lie in a single plane, whi ch is consistent 
with t he improved oblique pulsator model of iBigot fc Kurt3 
ll201ll) . KIC4768631, however, is unusual in having pulsa¬ 
tion maximum at maximum rotational brightness. Many 
roAp stars have pulsation maximum at rotational light min¬ 
imum, which can be associated with strong rare earth ele¬ 
ment spots. 

It tends to be the slow rotators, with periods longer 
than one month, that have aligned rotation and mag¬ 
netic axes llLmidgtreet^ MathvsilioOCll ) , or the more evolved 
stars, since North 1 19851) found essentially random angles of 
align ment for near-ZAMS stars (for a review see iMurphvI 
l2014h . The former is ruled out by the 5.2 d rotation pe¬ 
riod, but the latter case is a possibility since the non- 
adiabatic modelling suggests that the roAp-type pulsations 
would be expected for the star nearer the terminal-age main- 
sequence (TAMS) than the ZAMS. The luminosity, and 
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hence age, of KIC 4768731 is, however, uncertain dne to 
the lack of an independent distance determination. Further¬ 
more, Ap (SrCrEu) stars are supposed to take about half 
their main -sequence lifetime to develop their peculiarities 
||a^| 2009I1 . so the age is of interest. If Abt is right, then 
perhaps this star is not old enough to have accumulated 
very anomalous rare earth elements. Hence, KIC 4768731 
may prove to be an interesting test case for the develop¬ 
ment of chemical peculiarities in Ap stars. 
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